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Abstract 

New neutrino states Ub, sterile under the Standard Model interactions, can be coupled to 
baryons via the isoscalar vector currents that are much stronger than the Standard Model 
weak interactions. If some fraction of solar neutrinos oscillate into Uj, on their way to Earth, 
the coherently enhanced elastic //^-nucleus scattering can generate a strong signal in the dark 
matter detectors. For the interaction strength a few hundred times stronger than the weak 
force, the elastic ^-nucleus scattering via new baryonic currents may account for the existing 
anomalies in the direct detection dark matter experiments at low recoil. We point out that for 
solar neutrino energies the baryon-current-induced inelastic scattering is suppressed, so that 
the possible enhancement of new force is not in conflict with signals at dedicated neutrino 
detectors. We check this explicitly by calculating the iv induced deuteron breakup, and the 
excitation of 4.4 MeV 7-line in 12 C. Stronger-than-weak force coupled to baryonic current 
implies the existence of new abelian gauge group U(1)b with a relatively light gauge boson. 
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1. Introduction 



Standard Model (SM) of particles and fields must be augmented to include neutrino mass 
physics and perhaps extended even further to account for the "missing mass" of the Universe, 
or cold dark matter (DM). During the last decade the underground experiments [21 [3_1 HI 
El El [7J [8j [9] aimed at direct detection of DM in the form of weakly interacting massive 
particles (WIMPs) [1] have made significant inroads into the WIMP-nucleon cross section vs 
WIMP mass parameter space. Since no DM-induced nuclear recoil signal was found with the 
exception of two hints to be discussed below, they constrained many models of dark matter 
and ruled out some portion of the parameter space in the best motivated cases such as e.g. 
supersymmetric neutralino DM |10j . Higgs-portal singlet DM |TT] etc. 

It has been argued by some authors that although primarily designed to search for 
WIMPs, these experiments are in fact multi-purpose devices that can be also used for 
alternative signatures of other effects beyond Standard Model. In particular, using same 
instruments one can look for the absorption of keV-scale bosonic super- WIMPs [12] , search 
for the axion emission from the Sun [13], and also investigate some additional signatures of 
WIMP-atom scattering that exist in " non- minimal" WIMP models [2]. This paper extends 
this point further and opens a new direction: we show that neutrino physics beyond SM can 
also be probed with the dark matter experiments. 

Elastic scattering of neutrinos on nuclei pL5j is enhanced by the coherence factor N 2 , 
where N is the number of neutrons. Straightforward analysis [161 E] of the SM solar neu- 
trino elastic scattering rates on nuclei used in DM experiments reveal several basic points: 

• Despite the coherent enhancement, the scattering rates are way too small, leading to 
counting rates not exceeding 1CT 3 kg _1 day^ 1 keV^ 1 . Such low rates do not introduce 
any //-background to WIMP searches at the current levels of sensitivity. 

• The nuclear recoil spectrum is usually very soft, E r ~ (E u ) 2 / M Nud ~ few KeV or less. 

• Solar boron ( 8 B) neutrinos are the best candidates for producing an observable signal, 
because of the compromise between the relatively large flux and the energy spectrum 
extending to 15 MeV. 

Of course, at this point the DM experiments typically target much harder recoil and are 
far away from low counting rates induced by solar neutrinos. On the other hand, the last 
generation of dedicated solar neutrino experiments [HI UHl ED] have been extremely successful 
in detecting solar neutrinos via charged current reactions (CC), elastic scattering on electrons 
(ES), and Z-exchange mediated (NC) deuteron break-up [19]. It is the combination of all 
these three signals that led to a very credible resolution of the long-standing solar-neutrino 
deficit problem via the neutrino oscillation and the MSW mechanism [211 122] . 

However, it is easy to imagine that three active SM species z/ e , v T with their (almost 
completely) established mass/mixing parameters may not be the last word in the neutrino 
story. In this paper we consider a model of " quasi-sterile" neutrino z/& that has no charged 
currents with normal matter, no ES on electrons or other leptons, but has much enhanced 
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NC with baryons (NCB). We shall consider the strength of new NCB interaction to be much 
larger than Fermi constant, Gb ~ (10 2 — 10 3 ) x Gp. Such interactions can be mediated by 
new vector bosons of the gauged baryon number U(l)#, and for that reason we call this new 
hypothetical neutrino state as the "baryonic" neutrino v^. If any of the solar SM neutrino 
flavors oscillate into uj, within one astronomical unit, then the current DM experiments will in 
principle be able to pick it up via the coherently enhanced NCB signal. Whether such strong 
NCB would lead to a measurable energy deposition in the standard neutrino experiments 
requires special investigation and is addressed in this paper. We find that although the 
inelastic NCB scattering is enhanced by a huge factor G 2 B /G F , it is also suppressed by a tiny 
factor E^R%, where R n is a nuclear radius-related parameter. The resulting rate for NCB 
processes in neutrino detectors is then can be made comparable or smaller than the regular 
neutrino counting rates. 

It is somewhat tempting to relate the proposed z/f, model with the recently reported 
anomalies/signals in the direct DM detection. For a long time of course the DAMA and 
its successor DAMA/LIBRA experiments have been claiming jl] the annual modulation of 
the energy deposition in Nal crystals with the maximum in early June and minimum in 
December, which is consistent with the expected seasonal modulation of the WIMP-nucleus 
scattering rate. Then, last year the CoGeNT collaboration [5] has reported unexpected 
(in the null hypothesis) rise of their signal at recoil energies below E r =1 keVee. Given 
the mass of Ge nuclei, and typcial quenching factors in germanium, it is plausible that the 
rise of CoGeNT signal at low KeVee can be produced by resulting from oscillations of 
boron neutrinos, v$m ^b, and hypothesized enhancement of NCB can compensate for 
small neutrino scattering rates. It is also clear that mimicking DAMA signal with is also 
possible in a more restricted sense. Of course the usual seasonal modulation of the neutrino 
flux due to the eccentricity of the Earth orbit will have a minimum in the early July and a 
maximum in early January. However, the neutrino oscillation phenomenon is not monotonic 
in distance |23j, and if the oscillation length for vsm — > is comparable to the Sun-Earth 
distance, the annual modulation phase of the scattering signal can be reversed by n. We 
investigate this opportunity, and conclude that both CoGeNT and DAMA signal can be 
described with v^-type models (provided that DAMA data can tolerate a ~ 1 month phase 
shift). We further argue that if indeed this is the case, the model is very predictive, and 
there will be further ample opportunites for probing v\, both at DM and neutrino detectors, 
as well as at more conventional particle physics experiments. 

This paper is organized as follows: in the next section we introduce the class of i/& 
models and specify a parameter range that is the most perspective for the DM experiments. 
Section 3 addresses NCB elastic and inelastic scattering, including calculations of ^-induced 
activation of carbon and the deuteron breakup reactions. Section 4 studies the possibility of 
phase reversal in the annual modulation signal. Our conclusions are reached in Section 5. 

2. Baryonic neutrino and baryonic neutral currents 

The basic features of the model are as follows: we introduce a new gauge group U(l)#, and 
give all quark fields of the SM, Q, U,D, the same charge under U(1) B , which we call g&/3. 
We also introduce a new left-handed neutrino species v h that has a charge g\ under this 
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new group and no charge under any of the SM gauge groups. In the interest of anomaly 
cancellations it is also desirable to introduce a right-handed partner of v b with the same 
charge. Then the new group couples to the "vector-like" matter multiplets, and although 
SM+U(1)b will in general be anomalous, the anomalies can be cancelled at some heavy 
scales. Variants of this model may include some partial gauging of the SM lepton species 
under U(l)#. Neither the right-handed v b nor the details of anomaly cancellation will be 
important for this paper. Furthermore, we assume that U(l)# is spontaneously broken by 
the U(l)s-Higgs vacuum expectation value (4>b), and exactly how this happens will not be 
of direct consequence for us either. The relevant gauge part of the Lagrangian is then given 
by 

C = ~V%, + l -m\ V* + v bltl {id^ + gi V») Vb + J2 ^ II Psm + ^9blMq + C m . (1) 

i 

First two terms in ([T]) are the standard Maxwell-Proca terms for V^, the sum extends over 
all quark types and flavors, and D S m is the SM covariant derivative that includes gauge 
interactions appropriate for each quark species q. The mass part of the Lagrangian C m 
besides the usual SM mass terms should also account for neutrino masses, and generate 
mixing to a new state u b . In this paper we are not going to consider vector exchange with 
virtualities beyond O(10 MeV), and therefore it is convenient to switch from quarks to 
nucleons, 

gV^b 5Z<?7m<? ->■ 9bVp(jrytf) + n'ypn) + ... (2) 

Ellipses stands for 0(m]J~) terms associated with the V^Na^N part of the form factor, 
which will be small for any process we consider in this paper. The coupling of to the 
isoscalar vector current of nucleons = p^y^p + n^^n will have important implications for 
both the elastic and inelastic scattering of v b on nuclei. The exchange by the U(1)_b gauge 
boson creates the NCB Lagrangian, 

Cncb = Vblp.Vb f 9 * n Jj?\ (3) 
rriy + U p 

that in the limit of my 3> Q 2 is just a new contact dimension 6 operator with the effective 
coupling constant Gb- 

C NC B = G B xu bl ^ b J^; G B = 9 -^=Afx^-. (4) 

p m v GeV^ 

Here we have introduced an "enhancement" parameter M that quantifies how much stronger 
Gb is compared to the weak-scale value of 10~ 5 GeV~ 2 . We note that stronger-than-weak 
interactions among four neutrino species were considered earlier in e.g. Ref. [24J. The use of 
baryonic force as a mediator between SM and dark matter was considered recently in [25J. 

One may wonder if M as large as a 100 or a 1000 can be consistent with low energy data 
on meson decays, such as K — > tiv^v^. It turns out that due to the conservation of the baryon 
current, the loop amplitude for the underlying s — > dv b v h decay is additionally suppressed by 
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GpQ 2 , which compensates for all possible enhancements due to Af. (In contrast, the quark 
axial vector analogue of (DQ) will be strongly constrained to have Af < 1.) Thus, from quark 
flavor perspective, the baryonic portal (CQ) is one of the two "safe" portals (the other being the 
kinetic mixing with hypercharge [26]) that allow attaching stonger-than-weak interactions 
to the quark currents. We shall not pursue the meson decay constraints on the model any 
further in this paper, and turn our attention to the neutrino mass sector. 

The most natural way of having a UV-complete theory of neutrino masses is via the 
introduction of right-handed neutrinos states Nr. We can use the same singlet right-handed 
neutrinos coupled to the Higgs-lepton bilinears LH and fUggs^-neutrino z/& bilinears Vb^b i n 
a gauge-invariant way, 

C m = LHYN + v bL (j)hN + (h.c.) + ~N T M R N. (5) 

Here Mjj and Y are the familiar 3x3 right-handed neutrino mass matrix and Yukawa matrix, 
while b is the new Yukawa vector parametrizing the couplings of the left-handed part of Vb 
to N. Integrating out N states results in the low-energy 4x4 neutrino mass and mixing 
matrices, My, where i,j run over e,fi,r,b flavors. While of course a full four-state analysis 
can be done, we shall simplify our discussion by the following assumptions: 

1. The entries of 3 x 3 submatrix M act i ve , active will in general be somewhat larger than 
M active! & and Mb 5 b components so that the mixing pattern can be addressed sequentially: 
first the mixing of the SM neutrinos and then the admixture of the v^. 

2. A tri-bimaximal ansatz will be taken for the 3x3 mixing of the SM neutrino species 
for simplification, although having #13 = is not crucial. 

3. We shall assume a preferential mixing of Ub to ^2, with the relevant parameters that 
we call Ami an d @b, so that true mass eignestates are vi = cos Qb^i + sin ObW, vn = 
— sin QbVi + cos 6bVb- 

4. The sign of Gb will be chosen to ensure that the matter effects for Vb will not lead to 
the matter-induced z/ act i V c — > Vb transitions. 



The combination of these assumptions forms the following (simplified) picture of neutrino 
oscillations: inside the Sun the neutrino oscillations occur largely between u e and v + = 
(^ + v T )/y/2, 



Va- 




(6) 



while the "— " combination and 1% stay unexcited. We would need only the higher-end of 
the Boron neutrino spectrum where MSW effect dominates. Upon the neutrino exit from 

the dense region of the Sun, it represents an almost pure v 2 state, z/ 2 = \J\v e + \J\ v +i with 
individual flavor probabilities 

P e (Sun)~i; P + (Sun)~| P 6 (Sun) = 0. (7) 



4 



Then vacuum oscillations start building a non-zero probability for Vf, due to the vj and vjj 
being the true mass eigenstates in vacuum. Upon the arrival to the Earth, the following 
energy- dependent probabilities will approximate the neutrino flavor composition: 



P b (Earth) ~ sin 2 (2fl fe ) sin' 



P e (Earth) ~ - ( 1 - sin 2 (2# b ) sin 
P + (Earth) ~ M - sin 2 (2^) sin 



Am 2 b L(t) 
AE 
~Am 2 b L(t) 
_ AE 
Am 2 b L(t) 
AE 



(8) 



where L(t) is the Earth-Sun distance with a slight eccentricity modulation, 

~2Tr(t-t y 



L(t) ~ L Q ( 1 — ecos 



T 



L = 1.5 x 10 8 km; e ~ 0.167; t 



3 Jan. 



The most interesting range for Ami to consider is 



10~ 10 eV 2 < Am 2 b < A ™ 2 



Solar, atm' 



(9) 



(10) 



A scale of O(10~ 10 ) eV 2 is the so-called "just so" mass splitting that may introduce significant 
changes to the otherwise very predictable oc L~ 2 seasonal variations of the Vb flux at Earth's 
location. With Am 2 being much smaller than 10~ 5 eV 2 there is no danger of distorting 
KamLAND results [23J even for a relatively large angle 6b, although the matter effects for 
Vb could be significant. We also find it convenient to introduce the energy parameter E 
directly related to the mass splitting, 



E 



Am 2 b L Q 



An 



6.05 MeV x 



Ami 



lO" 10 eV 



2 ' 



ID 



that defines last zero of Pb as the function of energy, Pb{L = Lq, E = Eq) = 0. Since in all 



NCB rates Pb will enter in the combination with G 2 B , it is also convenient to define 



Af 2 



cS 



M 2 x - x sin 2 (20 b ), 



(12) 



so that in the limit of large Eq the oscillations average out and PbG\ ^ A/" c 2 ff x 10 10 GeV 4 . 

The pattern of masses and mixing considered here is not the most natural: we assume 
a pair of very degenerate vi and vu mass eigenstates replacing z/2 and Given that the 
mass of z/2, regardless of the hierarchy pattern, is always in between v\ and ^3, this will 
require some specific adjustments of the full 4x4 mass matrix. The search for more natural 
realizations of vsm ~^ v b oscillations with long oscillation length, including matter effects 
for a different sign of Gb, goes beyond the scope of this paper. The goal of the next two 
sections will be to find the sensitivity to A4fr in various processes involving the elastic and 
inelastic scattering of 

We would like to close this section with some model-building comments. A very intriguing 
question to ask is whether SM neutrinos would tolerate new large NCB. A conventional 
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answer is "no", as the so-called non-standard neutrino interactions (NSI) with quarks and 
charged leptons were addressed in a number of papers [27J and almost no room at 0(l)Gp 
level was found, let alone the much enhanced NCBs hypothesized in this paper. However, NSI 
studies [27] with rare exceptions [28] assume that the scale of the mediation is comparable 
to the weak scale, and ignore the possibility of light vector bosons communicating between 
neutrinos and baryons. As a matter of counterexample one could consider a model with two 
new gauge groups, U(1)b and the other being a quantized lepton flavor, e.g. or L T . The 
connection between two new vector sectors is given by the kinetic mixing term TjVffiV^ 



ft V 



Then additional effective interaction of a SM neutrino with the baryonic current is given by 

£eff USM1 " USM K V1 + D)« 2 + D) J * (13) 

Such interaction gives no contribution to the forward scattering amplitude and thus is not 
affecting neutrino oscillation, and it is 1/Q 2 suppressed in the large Q 2 regime, avoiding 
strong constraints from deep-inelastic neutrino scattering. It is then clear that the choice 
of my\, my 2 in the MeV range may allow having (Fl3l) at Q 2 ~ (1 — 10) MeV 2 to be 
considerably stronger than the SM weak force^. The interactions of type ( !T3|) can lead to the 
detectable recoil signal from elastic scattering of solar SM neutrinos, along the same lines 
as the ^-scattering idea advocated in this paper. A possibility of modified SM neutrino 
interactions such as (ITS]) can be very effectively tested using the proposed neutrino-nucleus 
elastic scattering detectors placed near the intense source of stopped pions [29J . 



3. Elastic and inelastic scattering of 

Elastic scattering of f& on nuclei will create a recoil signal regulated by the strength of NCB, 
and the probability of oscillation (jSJ). It can be picked up by the direct dark matter detection 
experiments with low recoil thresholds. Also, Vj, neutrinos can deposit a significant amount 
of energy on the order of a few MeV by activating excited nuclear states or via extra neutrons 
created by nuclear breakup. The main finding of this section can be summarized as follows: 
the ratio of the elastic to inelastic cross sections in the interesting neutrino energy range 
E v < 15 MeV, is governed by the following relation: 

^ fc - Nucl (elastic) ^ _A^_ ^ 
oVfc-Nuci (inelastic) E*R% 

where we took A ~ 100, Rjj 1 ~ 100 MeV, and E v ~ 10 MeV. It is this huge ratio that 
makes small-scale experiments such as [5] competitive in sensitivity to Uf, with the large- 
scale neutrino detectors. 

1 The author would like to acknowledge very stimulating discussions with B. Batell and I. Yavin on the 
possibility of the NSI enhancement. 
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3.1. Elastic scattering 



The differential cross section for the NCB elastic scattering of left-handed v\> on a nucleus of 
mass Mn with A nucleons is given by 

da el _ E 2 A 2 g 2 g 2 (l + cos6) ^ 1 x ^ + cos6) ^ 



d(cos6) 47r(M2+q2)2 ^ 

where the elastic scattering momentum transfer is q = (q 2 ) -1 ^ 2 = 2E sin(9/2) and cannot 
exceed twice the neutrino energy E. In the second relation we took My ^> E, which allows 
to shrink vector propagator. Using relations between the neutrino scattering angle 9, nuclear 
recoil energy E r , and the minimum neutrino energy required to produce i? r -recoiling nucleus, 



E r = — x(l-cos0); E"» = J*M>L, (16) 

771 jv V 2 



we rewrite the elastic cross section ( 115)) in the following form 



d ° Cl - 1 xG 2 B A 2 mJl-^^\. (17) 



dE r 2tt " V E 2 

One can readily see that the NCB cross section (fTTj) is related to the SM elastic neutrino- 
nucleus cross section [15] by G 2 B A 2 — > G 2 F (N/2) 2 substitution, where N is the number of 
neutrons (with small corrections in 1 — Asm 2 6w parameter). For the momentum transfers 
and nuclei considered in this paper the form factor corrections are < 5% percent and are 
ignored. 

Using cross section (fl7]) . the flux and the energy distribution of 8 B neutrinos [30] (hep 
solar neutrinos provide a small correction), we derive the counting rates as the function of 
interaction strength and the oscillation probability. For a moment, we neglect small seasonal 
modulations and take the limit of e — > . For the medium composed only of atoms with 
atomic number A, we approximate these rates by: 

^ * ^ x lsin 2 (2^)G|$s B x I(E r ,E ) 
ahi r 2ti 2 

recoils { A\ 3 Af 2 ff 

~85- — x — x-£xI(E r ,E ). (18) 



dayxkgxKeV \70 J 10 4 

The input (total flavor) flux of 8 B neutrinos is taken to be $sb = 5.7 x 10 6 cm~ 2 s _1 and 
ui jy oc Am p . 

The recoil integral I(E r , Eq) in eq. ([181) is given by the convolution of the 8 B energy 
distribution, the energy-dependent part of oscillation probability, and the kinematic factor 
in the cross section reflecting neutrino helicity conservation: 



poo / /77imin\2\ 

I(E r ,E ) = dE 1 - ^— J- x f SB (E) x 2 sin 2 

JE™ in (E r ) V & J 



ttEo 
E 



(19) 



Here the distribution function is normalized as J u E fs B (E)dE = 1. For the limit of large Eq 
(fast oscillations), the last multiplier in (fT9l) becomes 1. If a detector threshold corresponds 
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Figure 1: Expected recoil event rate in Germanium in units of recoils/day/kg/keVee as the 
function of E r in keVee. The NCB enhancement factor, A4ff = 100. A, B and C lines 
correspond to E = oo, E = 12 MeV and E = 14 MeV. 

to recoil energies that require E mm to be above the end-point of 8 B neutrino spectrum, 1 = 
(apart from small corrections from hep and diffuse supenova neutrinos). It is the case for 
most of the existing WIMP detectors, but not for all of them. 

In real detectors registering ionization such as [U |5] it is the electron equivalent of the 
energy release rather than the recoil energy that is detected. We take the relation between 
the two by following recent DM-related analyses [3T1 132] 



The second relation is far less precise than the first one [32J. 

The counting rates in germanium resulting from scattering of Ub created by the oscillations 
of 8 B and hep solar neutrinos are presented in Fig. 1. We have taken three cases of mass 
splitting: E > E™ x , and E = 12, 14 MeV. The NCB rates are plotted for the value of 
N"£ s = 10 4 . For this enhancement factor, the resulting counting rates are clearly within reach 
of current generation of low-threshold germanium detectors (i.e. CoGeNT). 

Inspection of Figure 1 shows that choice of different mass splitting that makes the oscil- 
lation length comparable to 1 A.U. influences the shape of the spectrum. This is because the 
most important part of the spectrum for the recoil in excess of 0.5 keVee is above neutrino 
energies of 10 MeV, where the 8 B neutrino spectrum is already sharply falling. The neutrino 
oscillations with Eq close to 12 MeV will lead to the suppression of higher E r and to the steep 
rise of the signal at lower E r . The sharp end of the neutrino spectrum prevents other Ge 



Ge : £ r (keVee) ~ 0.2 x (£ r (keV)) L12 
Na in Nal : £ r (keVee) « 0.33 x E r (keV). 



(20) 
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experiments with higher threshold like CDMS [2] to probe the NCB scattering in the regime 
of large recoil where [2] has strong sensitivity. The signal from the recoil due to v b neutrinos 
is very similar in morphology to that of sub-10 GeV scale WIMPs. This is because a typical 
momentum transfer in a heavy nucleus - light WIMP collision is q ~ m W i mp v ~ 10 MeV, 
which is about the same for 8 B neutrino scattering. There is one kinematic difference though: 
the back-scattering of WIMPs that produces hardest recoil is kinematically allowed, while 
for neutrinos it is forbidden by helicity conservation. This additionally limits the capabilities 
of high-threshold experiments to detect v b neutrinos in comparison with light WIMPs. 

Is it possible to use v b as another speculative explanation of CoGeNT results [S]? The 
overall event rate can indeed be reproduced well with Af c{i ~ 10 2 — few x 10 2 , depending on 
Eq. For the large Eq parameter, the enhancement factor of M e s = 10 2 seems sufficient: it 
gives 7 recoils/day /keVee at E r = 0.7 keVee, which is about the same as the experimental 
data suggest after accounting for the efficiency [51 [32] . The shape of the predicted signal is 
also similar to the counting rate profiles observed by CDMS at the Stanford Underground 
facility [33]. Fitting the exact spectral shape of excess events at CoGeNT falls outside the 
scope of our current investigation. We should also note that the expected total counting 
rate for the material used in CRESST detectors [6] due to the neutrino-oxygen scattering is 
given by 

recoils A/"^ 

R in CaW0 4 (£r > 10 keV) ~ 0.2 X — X -f , (21) 

day x kg 10 4 

which is well within their detection capabilities for A/" c ff ~ 100. Other methods in devel- 
opment that use liqud helium as a detecting medium with a potentially very low energy 
threshold [31] also look promising for detecting ^-induced recoil. It is also important that 
the choice of very low-mass target such as 4 He will allow discriminating between > 5 GeV 
WIMPs and i/&'s: the effective recoil energy goes down at Mn < Mwimp, while it becomes 
larger for v b scattering. 



3.2. Inelastic scattering 



Unlike light WIMPs that can cary significant momentum but very little energy, v b can easily 
lead to an MeV-scale energy deposition. Here we turn our attention to the NCB inelastic 
processes and will adress the following issues: the NCB deuteron breakup, and the NCB 
excitation of the first 2 + resonance in 12 C resulting in 4.4 MeV 7 line: 

d + u b u b + n + p (22) 
12 C + u b v b + 12 C*(4.44 MeV) — )■ v h + 12 C + 7 (23) 

The main scientific question to answer is whether the enhanced values of G\P b can be 
consistent with the constraints provided by SNO on "extra neutrons" from ( 1221) and by 
Borexino and other liquid scintillator detectors on "extra gammas" from (T23T) . There are of 
course other processes that one has to consider in a more comprehensive study, including the 
excitation of oxygen, the breakup of 13 C to 12 C+n etc, but they will all follow the scaling 
in Eq. ( fl4l) . The earlier studies of the nuclear excitations due to the differenty type of 
neutrino couplings can be found in [35]. There are also elastic channels of energy deposition 
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via u b + p — > v^ + p [36J , but the proton recoil from the scattering of 8 B neutrinos would fall 
below the detector thresholds. 

To understand the origin of ratio ( TT4l) one does not have to perform any sophisticated 
calculations. We consider the scattering of ~ 10 MeV energy neutrinos, so that their wave- 
lengths are much larger than the characteristic nuclear size of a few fm. Therefore, one can 
safely expand the nuclear matrix elements in series in q, or in neutrino energy E, as q is 
bounded by E. Here is how the inelastic matrix element of the /i = component of the 
isoscalar vector current between the deuteron bound state and np continuum will look 
like in this expansion: 

(d\ exp(iqr (n) ) + exp(iqr (p) )|np) (24) 
= 2(d\np) + iq ■ (d\r^ + r^\np) — -^-{dlr^r^ + r^r^\np) = —^-^{d\r^ri\np)^ 

where r", r^' are the position operators for the neutron and the proton. The zeroth and 
first order terms in q = |q| are trivially zero due to the orthogonality of the wave functions 
( r («0 _|_ r (p) i s the center-of-mass operator and cannot mediate inelastic transitions). In the 
last line we have introduced the relative position vector r = r^ n ^ — r^ p \ and the quadratic 
in r operator can be further separated into the isotropic "charge-radius" and quadrupole 
components. For the + — > 2 + transition in 12 C only the quadrupole part will matter. It 
is of course instructive to revisit the SM deuteron breakup [37], and observe that isoscalar 
vector component of the standard weak current gives a very minor contribution to the total 
cross section at low E u due to this q 2 suppression of the amplitude. The SM rate is of course 
dominated by the isovector axial- vector current that corresponds to the difference of nucleon 
spins s^") — s^, an operator that has non-zero inelastic matrix elements even in the 0(q°) 
order. 

We perform the calculation of the NCB-induced deuteron breakup using the "zero-radius" 
approximation of the initial and final state wave functions, 

^in(r) = exp(-M-); ^f( r ) = exp(ipr) (25) 

V47rr 

K = y^2E b p: = 45 MeV; p 2 = p 2 = 2fi(E - E f - E b ). (26) 

Here E and Ef are the initial and final energy of z/&, E b = 2.2 MeV is the absolute value of the 
deuteron binding energy, and fi ~ (m n + m p ) /4 is the reduced mass of the proton- neutron 
two-body system. Relative momentum p of the final state is fully determined from the 
neutrino energies, as the recoil of the deuteron center-of-mass is negligible. Parameter k is 
the familiar "bound state momentum", related to the inverse size of the deuteron, k ~ R^ 1 , 
and its relative smallness reflects large spatial extent of the deuteron. In a slightly excessive 
for the problem at hand language our calculations correspond to the leading order of the 
pionless effective field theory EH]. They can be systematically improved if needed, or 
treated with the more elaborate nuclear physics tools (see e.g. [39J). None of this will of 
course change the order of q in which the effect first occurs. 

Straightforward calculations give the differential over the final neutrino energy cross 
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section: 

d&d—tnp 



G 2 E 2 m r , 



K 5 p 



dE 



f 



87T 2 (p 2 + K 2 



5k 4 



2 „ 10 

3 EE ' + T 



2 
3 J 



£' J £'7 - -BE} + Ej) 



(27) 



The result for dad^ np shows a 0(E a k~ 4 ) suppression in agreement with 0(q 2 ) of the deuteron 
matrix element and in agreement with (Q3J). Judging by the size of the subleading corrections 
in the SM calculations [37] we expect this answer to hold within ~ 20% accuracy. 

The final integral over Ef in the interval from to E — Ej, gives the total NCB deuteron 
breakup cross section. In Figure 2, upper pannel, we plot cr d ^ np for the Standard Model 
neutrino and for Vb with the choice of enhancement factor M 2 = 10 4 . As expected, the NCB 
cross section has a faster rise with neutrino energy due to quadrupolar nature of the NCB 
interaction. In Figure 2, lower pannel, we also show the convolution of the cross section with 
the energy distribution of boron neutrinos times the energy-dependent part of the oscillation 
probability, 2 sin 2 (7ri?o/-E'). As in the previous subsection, the NCB rates are considered for 
large Eq and for Eq = 12, 14 MeV, while the enhancement factor is kept at Af c s = 100. 
The areas under curves give the total effective cross sections, and when multiplied by $s B , 
correspond to the breakup rate per deuteron. For the three NCB cases considered here have 
the following comaprison to the SM rate: 



oncb Af 2 

r^j 

CSM NC 10 4 



cff 



x (0.14, 0.06, 0.13) at E = oo; 12 MeV; 14 MeV. 



(28) 



A 15% increase in the neutron production rate at SNO can be tolerated, and if one chooses 
a sizeable 9b so that the active neutrino flux is slightly less than SM+SSM predict, the total 
neutral current rate may not even change. We conclude that SNO NC events leave enough 
room for the possible Af 2 s ~ O(10 4 ) (and slightly higher) enhancement of the NCB rate. 

We now address the 12 C — > 12 C* activation due to Vb- To avoid the complications arising 
from nuclear physics we shall assume that both the ground state and the first excited state 
of 12 C are given by 3a configurations. This is a very well justified assumption, which leads 
to a relation between the matrix elements of the baryonic current and electric current, 

(E 2 -t 



(0 + |J i (0) |2 + ) = 2(0+|Jr|2 + ) = 2- 



6 



Eo+ k( +|Q-|2+> 



(29) 



Only the lowest order in q terms are retained here, and the fi = component can be 
restored from gauge invariance. The factor of 2 in (T29]) comes from the fact that the baryonic 
charge of a-particles is twice larger than its electric charge. The information on the value 
of the transitional quadrupole moment (Qij) can be extracted from the 12 C* decay width 
r = 1.08 x 10" 2 eV: 

qnr 

Kg ^ |2= aTAEF = (3 - 3fm)4 ' (30) 
where the value of quadrupole moment squared is averaged over the arbitrary projection of 
the 2 + angular momentum, and AE = E 2 + — E + = E — Ef = 4.439 MeV. Note that we 
define electric quadrupole and electric current without the " e" , and account for a explicitly 
in (1301) . The total inelastic cross section for the ^-induced 12 C — > 12 C* transition is given by 



<Tl2 C _>.12 C * 



8G 2 B E*{E - AE)\(Q?p 
81tt 



39 o 19 o 39 a 9 r 

ix H x x H x x 

8 4 16 16 



(31) 
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Figure 2: Top pannel: Deuteron breakup cross section for the SM NC processes (top curve) 
and for the NCB ^-neutrinos (bottom curve). The NCB cross section is plotted for N 2 
enhancement factor of 10 4 . Bottom pannel: same cross sections convoluted with 8 B energy 
distribution and the energy dependent part of the oscillation probability. The top curve is 
the SM NC distribution of effective cross section, and A, B, C are the same for the NCB 
with Eq = oo, 12, 14 MeV and M 2 S = 10 4 . The areas under the curves give the proportion 
of neutrons produced via SM NC and NCB processes. 
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where x = AE/E. The benchmark value for this cross section at E = 8 MeV and M = 1 is 
2.5 x 10~ 48 cm 2 . 

With this cross section the effective rate of injection of 4.4 MeV gamma quanta in pseu- 
documene (scintillating material used by the Borexino experiment) is estimated to be 

R(4A MeV) ~ (0.05 - 0.15) x 7 injeCti ° n ; x % (32) 
v J K ' 100 tons x day 10 4 v ; 

This is not a large rate by any measure, but it is nevertherless comparable to the counting 
rates in 3-5 MeV window from 8 B ES processes and from the 208 T1 background events [40J. 
The actual counting rate should be obtained by applying to (|32|) the efficiency factor that 
the collaboration can extract from their calibration data and simulations. At this point we 
can only conclude that there must be some sensitivity to NCB at A4s ~ 10 2 level at the 
large scale neutrino detectors that use carbon-based scintillators. More definitive statements 
and perhaps stronger sensitivity to Ub can be derived from dedicated analyses. Moreover, the 
search for the extra 7-lines in a different energy range was already performed by the Borexino 
collaboration in connection with hypothetical Pauli-forbidden 12 C decays [41J. The search for 
NCB would represent a far less exotic physics cause in our opinion. One could also conduct 
similar searches of ^-induced excitation of 16 O nuclei using SNO and SuperK data. 



4. Annual modulation of z/& rates 

In this section we would like to address the question of seasonal modulation of the NCB rate. 
The seasonal modulation of solar neutrino rate was observed by SNO and SuperK collabo- 
rations [HI H2]. It exhibits full agreement with expected oc L~ 2 (t), 3.3% modulation of the 
neutrino flux, with an appropriate minimum in the summer (Northern hemisphere). The hy- 
pothetical NCB elastic scattering rate will have the same modulation pattern as long as Am 2 
is large, or in other words at E ^> E v so i ar . In the opposite limit of E <C E v so i ar the mod- 
ulation effects are suppressed because oc L^L^it) sin 2 [nEoL^^LoE)^ 1 ] — > (tiE q /E) 2 , 
which is time- independent. 

However, it is easy to imagine that the flux of neutrinos can have a more intricate 
seasonal modulation pattern. For example, if E is between the maximum of the 8 B neutrino 
spectrum and its end-point, the high-energy fraction of the distribution will have a higher 
flux in the summer. This is best illustrated in Figure 3, where the expected flux of i/& 
resulting from oscillations of boron neutrinos is convoluted with the time- dependent part 
of P b , LlL~ 2 {t) sin 2 ''[nEoL^iLoE)- 1 ] at E = 12 MeV. The two curves correspond to t = 
^perihelion (~ 3 Jan) and tapheiion (~ 4 Jul). Although on average there is more v\, neutrinos 
arriving to the Earth in January, in the most relevant range of energies, E > 10 MeV, the 
flux in July is larger. Therefore, for this fraction of neutrinos there is a phase reversal, and 
the elastic scattering rates will reflect that. 

Next we calculate the expected seasonal modulation in the counting rate, 

dE r ~ 2 \dE~ r 



Jul 



dR 
IK 



Jan 



(33) 
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Figure 3: Boron v\> neutrino flux modified by the time-dependent part of the oscillation 
probability 2L 2 1 L" 2 (t) ^[kEqL^LqE)- 1 } with E = 12 MeV. The black curve is for July, 
and the gray curve is for January Although the total integral under the gray curve is bigger 
than under the black one, it is the high-end of the spectrum that would determine rates at 
the existing DM detectors, where July rates are larger. 

for Nal detectors using the quenching factor from Eq. (1201 . We would like to remark in 
passing that for some ranges of neutrino energies there can be a significant departure from a 
simple time-sinusoidal function, but to observe such effects one would probably require very 
high statistics and very good energy resolution. Modulation rates, dR mo d/dE r , as defined in 
Eq. ( 133|) are plotted in Figure 4 for the same three choices of E and Af c s = 100 as before. 
One can see that indeed modulation of both signs is possible, and that the rate of modulated 
NCB signal at this A/" e fr is indeed probed by DAMA/LIBRA experiment [4], which is sensitive 
to modulation amplitudes O(10 _2 )cpd/kg/keVee. 

Is it possible that z/^-Na elastic scattering is behind the DAMA/LIBRA seasonal modula- 
tion anomaly? The magnitude of the predicted modulation can be in a very good agreement 
with DAMA results [4J. Moerover, as we saw in the previous section, A/" 2 ff ~ 10 4 is thus 
far consistent with other observations and constraints (and with simultaneous explanation 
of CoGeNT low E r anomaly). Of course the phase of the DAMA results will require E Q to 
be in the right range. But even then, would the early July maximum be consistent with 
DAMA/LIBRA claims of the oscillation phase? The best fit point for the maximum is about 
4-5 weeks different from the t ap h e iion 0]- It would be interesting to know if the early July 
maximum is actually excluded by DAMA data, and the criticism expressed in Ref. [13] 
about the errors on the phase being too tight be properly addressed by the collaboration. 

Going away from idea, one can also notice that many other exotic physics explanations 
of DAMA signal can be invoked (if it allows to tolerate ~ 1 month phase shift). For example, 
the emission of solar axions with their subsequent absorption in DM detection experiments 
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Figure 4: Modulation of the counting rate in recoils/kg(NaI)/keVee for ^-scattering on Na. 
As before, A curve is for large Eq, B is for Eq = 12 MeV and C is for Eq = 14 MeV while 
A/" c 2 ff = 10 4 . Both signs of modulation are possible. 

can be a cause of low-energy ionization signal [13]. On this picture, one can super-impose 
the oscillations of axions into some "sterile axions" with the oscillation length similar to 
Lq in order to break the monotonic L(t)-dependence, and flip the phase of the modulation, 
achieving results similar to those of Figure 3 and 4. 

5. Discussion and conclusions 

The oscillation of SM neutrinos to some new neutrino state on the way from the Sun to the 
Earth is a realistic possibility. We have shown that there exists the whole class of models 
where neutrino physics beyond SM can be probed by the low-threshold DM detectors, which 
become equally sensitive or even more sensitive to this type of neutrinos than the large 
scale neutrino detectors. Such models require that new neutrino states v h (or modified 
SM model neutrino interactions in the spirit of Eq. f fl3|) ) couple almost exclusively to 
the baryon current. The isoscalar vector properties of this current lead to a very strong 
enhancement of the elastic over inelastic scattering, <T e i ast i C / ^ inelastic ~ 10 8 , providing an 
unexpected competitiveness factor to small scale experiments such as CoGeNT. We have 
shown the effective strength of the NCB can be much larger than the weak-scale value 
without being in conflict with any of the observational data. 

We have also shown that the recent anomalies in direct DM detection, such as DAMA 
and CoGeNT, can be explained by the z^sm — > oscillation of 8 B neutrinos with subsequent 
scattering of v h on Ge and Na nuclei. (This statement relies on the assumption of £ a pheiion 
being consistent with DAMA/LIBRA modulation phase.) This may look counter-intuitive 
at first, but we have shown that the phase flip of seasonal modulations is possible for the 
high-energy end of the 8 B spectrum. This is a very speculative explanation (and perhaps 
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equally speculative as the WIMP recoil explanation), but it is interesting enough to motivate 
further studies. In particular, we believe that the Borexino collaboration can perform the 
search of the //^-activated 4.4 MeV line in 12 C, and probably surpass the sensitivity to the 
NCB enhancement factor Af c g of 100. At the same time it seems apparent that further 
technological developments of low-threshold WIMP/V& detectors are required. Should the 
current low-energy anomalies in DM detectors firm up to constitute a definitive signal of new- 
physics-induced recoil, some significant efforts and different mass targets would be required 
to observationally distinguish between the low-mass WIMP and signals. 

Below, we would like to discuss further implications of the models involving new neutrino 
states with enhanced baryonic currents. 

• Collider implications. If the Gb ?t 100GV, and if the new interaction is truly contact, 
the proton-antiproton collisons will lead to strong new sources of missing energy signals 
in f&z?& channel and will most likely be excluded by the Tevatron experiments. This will 
not happen, however, in models of relatively weakly coupled mediators with sub-GeV 
mass. Therefore the collider searches should be able to place an upper bound on my. 

• Fixed target implications. A GeV-scale U(l)# baryonic vector, the carrier of NCB 
interaction, can be produced in the collisions of energetic proton beams with a target. 
Immediate decays of these vectors will generate a flux of state that can be searched 
for at near detectors via their NCB interactions. This is very similar to the ideas of the 
"MeV-scale DM beams" discussed previously in [44]. There can be also implications 
for the terrestrial anti-neutrino physics, as matter effects induced by ^-exchange can 
be large for antineutrinos [2S] . Enhanced neutral currents of i/& neutrinos may help 
explaining the long-standing puzzle of the LSND anomaly [45], perhaps borrowing 
some elements of the recent suggestion [46]. It also has to be said that over the 
last two years a lot of efforts have been invested in systematically searching for the 
"kinetic mixing" (or hypercharge) portal (see e.g. (47] and references therein). Barionic 
portal is another example of a perfectly safe from the model-building perspective way 
of introducing stronger-than-weak forces at low energy, and therefore it should be 
systematically searched for using proton-on-target facilities. But perhaps the most 
NCB-search effective type of experiments to perform with proton beams is the proposal 
[29] of a neutrino-nucleus elastic scattering detector. 

• Cosmological implications. A new light neutrino state (and two neutrino states if 
right-handed copy of Vb is also light) can be at borderline of what is allowed by early 
cosmology and observations of light elemental abundances (most recent analysis can be 
found in |48j). Is the model with much enhanced baryonic currents has a chance to be 
consistent with these constraints? Actually, despite the interaction strength of 100 or 
1000 of Gf-i Vb will decouple from thermal plasma earlier than the SM neutrinos. That 
is because its thermalization rate will be proprtional to the baryon-to-photon ratio, 
which is a small number O(10 -10 ). Therefore actual decoupling of Vb may happen with 
the decays and annihilations of abundant hadronic species at temperatures of ~ 100 
MeV, and therefore BBN bounds from over-population of radiative degrees of freedom 
can be easily evaded. 



16 



• Astrophysical implications. Another interesting aspect of v\> models is their NCB pro- 
duction in stars. In the SN Uf, will not provide new effective energy sinks because they 
would not escape freely the explosion zone. However one should expect that a com- 
parable to the SM number of f& neutrinos is created, so that one could detect them 
using the same DM/Vb detectors. Should a nearby SN explosion happen, the existing 
neutrino scintillator detectors can pick up the z/^-NCB signal that would appear as a 
much enhanced v^, v T -NC signal considered in [49] (modulo the uncertainty in effective 
temperature for u b ). 

• Rare decay implications. Relatively large NCB currents should open new channels 
for the missing energy decays of B and K mesons. As argued in this paper, the 
conservation of the baryonic current makes it a relatively safe portal compared to e.g. 
scalar or axial-vector portal. Nevertheless, if the K — > ttV decay is kinematically 
allowed, it may lead to the underlying two-body signature of K — > n plus missing 
energy decays, making it an appealing target for the next generation of precision kaon 
physics experiments. 
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